Electronic structure of fluorides: general trends for ground and excited state properties 
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The electronic structure of fluorite crystals are studied by means of density functional theory within the local 
density approximation for the exchange correlation energy. The ground-state electronic properties, which have 
been calculated for the cubic structures CaF2,SrF2, BaFi, CdF2, HgFi, p-PbF2, using a plane waves expansion 
of the wave functions, show good comparison with existing experimental data and previous theoretical results. 
The electronic density of states at the gap region for all the compounds and their energy-band structure have been 
calculated and compared with the existing data in the literature. General trends for the ground-state parameters, 
the electronic energy-bands and transition energies for all the fluorides considered are given and discussed in 
details. Moreover, for the first time results for HgF2 have been presented. 

PACS numbers: 71.20.Ps, 71.15.Mb, 78.20.Bh, 71.15.Qe 



0^ 

(N 






a 



I 

o 



> 

o 



X 



I. INTRODUCTION 

Fluorides and fluorite-type crystals have attracted much in- 
terest for their intrinsic optical properties and their potential 
applications in optoelectronic devices. The technological im- 
portance of CaF2 is due mainly to its optical properties; in- 
deed CaF2 has a direct band gap at F of 12.1 eV and an in- 
direct band gap estimated around 11.8 eV. ' Calcium fluorite 
CaF2 as well as all the other fluorides, is a highly ionic insu- 
lator with a large band gap, and its lattice structure is a cubic 
Fniim within three ions per unit cell, i.e. one cation placed 
in the origin and two anions F are situated at±(|fl, |fl,|fl)." 
Here we shall consider CaF2, SrF2, BaF2 (with cations be- 
longing to the II group) and CdF2, HgF2 (with cation belong- 
ing to group IIB) and, finally, jS-P/jF? (with cation belonging 
IV group). We refer to them respectivelly in the present text 
as Il-compounds, IIB-compounds and ^-PbF2. We therefore 
propose the first systematic study of the electronic properties 
of some fluorides compounds within the same computational 
approach. Until now the theoretical studies of these com- 




(a) conventional cell 



(b) Brillouin zone for Fm3m 



Figure 1 : (a) Fluorite structure. Small spheres represent F' 

ions and the larger ones the cations++. (b) The Brillouin 

zone of space group Fniim. 



pounds have been tackled in separate forms and with different 
techniques. Therefore no general and systematic trend for the 
family of these compounds could be obtained. The DFT-LDA 
studies are of particular importance and are benchmarks for 
subsequent researches to perform excited states and optical 
properties calculations. The latter rely on more sophisticated 
techniques which must start from well converged DFT-LDA 
calculations (e.g. perturbative GoWg, self-consistent GW, BSE 
etc.). ' ' In the following lines we proceed to consider the ex- 
perimental data and the theoretical results for fluorides present 
in the Uterature. 

a. Experimental Scenario Due to its importance in ap- 
plication and basic research, experiments on CaF2 and flu- 
orides compounds have been carried out for a long time. 
The optical constant of CaF2 in the extreme ultraviolet 
has been studied by discharge-tube technique.' Reflectance 
spectra, transition energies for CaF2, BaFo, SrF2 and other 
ionic compounds have been measured by synchrotron radi- 
ation facility later on.' Studies on l5-PbF2, CaF2, SrF2 e 
BaF2 dielectric properties have been published as a func- 
tion of pressure and temperature by capacitance and dielec- 
tric loss measurements. Satellites in the X-ray spectra for 
CaF2, SrF2, BaF2 , density of states for intraband transitions, 
have been studied by photoelectron spectrometry.^ The ef- 
fects of Eu defects in P-PbF2 {PbF2 : Eir'^) relative to fluo- 
rescence/electronic excitation spectra and dielectric relaxation 
have been analyzed by laser absorption spectra techniques.^ 
j5-PbF2 and CdF2 mixed crystals absorption coefficients have 
been reported by spectrophotometry measurements. Neutron 
diffraction techniques have been employed to determine the 
j3 -PbF2 pressure and temperature phase diagram. ' a — PbF2, 
j3-PbF2 and others compounds absorption spectra and elec- 
tronic transitions have been studied by polarized reflectivity. ' ' 
Polarized deep and vacuum ultraviolet light measurements 
permitted the study of birefringence of CaF2 and BaF2. The 
study by different techniques of the phase diagram of HgF2 
and other Hg compounds should also be mentioned here'' 
(even if for HgF2 no results appear in this work due to hydra- 



Table I: Optimized lattice constants of fluorides. In columns 

"LDA" and "Theory" we show previuos theoretical results 

(DFT-LDA and others respectively), while in column "Exp." 

data after different experimental references are reported. 



ao[A] 


present 


LDA 


Theory 


Exp. 


CaF2 


5.30 


5.34"' 


5.46' ' 


5.46*' 


SrF2 


5.68 


- 


5.79'^ 


5.78'^^ 


BaF2 


6.09 


- 


6.26^2 


6.1715 


CdF2 


5.31 


5.36 


5.39'' 


5.46^" 


HgF2 


5.47 


- 


5.55--' 


5.54^^ 


li-PbF2 


5.77 


- 


5.94 


5.94 



Table II: Bulk moduli of fluorides. Data in column "Theory" 

are previous theoretical results. Data in column "Exp." are 

after different experimental references. 



So[GPa] 


present 


Theory 


Exp. 


CaF2 


103.01 


91 


90-82 ■■ 


SrF2 


83.75 


- 


- 


BaF2 


69.39 


50^^ 


5946 


CdF2 


123.96 


123-"- 


- 


HgF2 


117.03 


- 


- 


l3-PbF2 


93.22 


60 


64 



II. COMPUTATIONAL DETAILS 



tion of the sample). Recently schematic experimental phase 
diagrams for HgF2 and HgF have been also reported. '^ 

b. Theoretical Scenario Various theoretical methods 
have been applied to study either the ground state or the ex- 
cited states of the fluorides compounds. One of the first works 
of relevance is the one in which elastic constants, pressure 
derivatives of 2'"' order elastic constants, static dielectric con- 
stant and its strain dependence for CaF2, SrF2, BaF2, have 
been calculated within a shell model. '^ The energy bands 
and reflectance spectra of CaF2 and CdF2 have been deter- 
mined afterwards within a combined tight-binding and pseu- 
dopotential method. ' Mixed crystals ofCaF2, SrF2, CdF2, (5- 
PbF2 have been studied with respect to their energy bands and 
DOS within LMTO method. ' A phenomenological method 
has been then applied to calculate specific heat of j5-PbF2, 
CaF2 and other compounds.'^ Linear and non-linear optical 
properties of the cubic insulators CaF2, SrF2, CdF2, BaF2 
and other compounds have been determined by first princi- 
ples OLCAO. Points defects study in CdF2 have been per- 
formed within the plane wave pseudopotential method. - With 
respect to electronic excitation properties and energy band- 
gaps, electronic band structure of CaF2 and other compounds 
have been determined by DFT-GW, using plane wave basis 
set and ionic pseudopotentials (PW-PP) scheme."' Using the 
hybrid B3PW functional, the electronic structures of defected 
fluorides, namely CaF2 and BaF2, have been calculated."""""' 
The £2(0) function after an iterative procedure using an ef- 
fective Hamiltonian has been calculated for CaF2 and GaN,-^ 
within PW-PP considering a screened interaction for e-h cou- 
pling. Native and rare-earth doped defects complexes in j3- 
PbF2 have been studied by atomistic simulation."" 



All the calculations for the fluorides under study have been 
performed using density funcional theory (DFT) method 
implemented in the plane-wave basis code VASP"^'-' Projec- 
tor augmented wave pseudopotentials (PAW)""'"' have been 
used in the localized density approximation (LDA) for the ex- 
change correlation energy treated within the scheme of Ceper- 
ley and Alder parametrized by Perdew and Zunger.^'-'^^^^ 

Relativistic effects have been included in the 
calculations-'"*'-''' ""''"'' and spin-orbit coupling has been 
considered for the valence electrons. 

The conventional cell of the crystals is shown in Fig. la, 
in which the ions F'^' (drawn as little spheres) form a cu- 
bic sublattice surrounded by a faced cubic center lattice of 
cations*++' (in the picture shown as large spheres labeled e.g. 
as Calcium). Fluorides with cations belonging to the II and 
IIB groups show a stable phase in this crystallographic struc- 
ture. On the other hand, at low pressure the PbF2 show two 
structural phases, namely orthorhombic (a) and cubic (j3). 
Although the cubic phase l5-PbF2 is the most stable in am- 
bient condictions, while the orthorhombic a-PbF2 becomes 
stable at high pressure. ' 

In Fig. lb, the Brillouin zone has been drawn for the Fniim 
space group; the paths in the fc-space chosen in our calcula- 
tions for the electronic band structures have also been shown. 
All calculations are performed with the total energy conver- 
gence within 1.5 • lO^^eV with kinetic energy cut-off depend- 
ing on the cation of the compound under study (at least of 
550ey). A Monkhort-Park fc-point mesh of 4 x 4 x 4 has 
been chosen to relax the cell parameters, till the largest value 
for the interatomic forces result smaller than 1.5 • 10^^ eV/A. 
Lattice constants for each fluorite (ao in Tab. I), as well as 
the bulk moduli (Bo in Tab. II) have been calculated with the 
Murnaghan equation of state. "'^ 



In this paper, we are interested in either the structural and 
the electronic properties of each fluorite, and in compari- 
son/trend studies for the whole crystallographic family. We 
have computed therefore the electronic and structural prop- 
erties of different fluorides, CaF2, SrF2, BaF2 (with cation 
belonging to group II), CdF2, HgF2 (with cation belonging 
to group IIB), j5-PbF2, within the same first-principles pseu- 
dopotential plane-wave method. 



III. RESULTS AND DISCUSSION 

Trends for lattice constants and bulk moduli are showed in 
Fig. 2 with the corresponding experimental data from litera- 
ture, if available. From Table I and Fig. 2 an overall good 
comparison appears between experimental data and present 
results relative to lattice constants with maximum deviation 
of 3%. Compounds with cations belonging to the II and IIB 



Table III: For the fluorides considered here, bandwidths and energy gaps as found in our DFT-LDA study. Direct and indirect 
DFT-LDA band gap energies. In the first column, current results are shown. Data are in e V. 





Minimum Direct Band gap 


) 


Minimum Indirect Band gap 


Valence Bandwidth 


Solid 


present LDA GW 


Exp. 




present 


LDA GW Exp. 


present LDA GW Exp. 


CaF2 


r^r 7.71 7.11 11.80 


12.10 


x^r 


7.43 


6.85'^ 11.50-^** 11.80' 


3.17 2.82^** 3.26-"^ 3.20 ' 


SrF2 


r^r 6.99 


11.25' 


x^r 


6.89 


6.77''^ 10.60' 


2.33 2.80' 


BaF2 


r^r 6.67 


11.00' 


{llo)^r 


6.58 


7.19'** 10.00 


1.78 2.50' 


CdF2 


r^r 3.37 3.30-" 




w ^r 


2.94 


2.80-" 7.80-''-' 


5.79 


HgF2 


r^r 0.35 




r ^L 


4.16 




6.38 


P-PbF2 


X^X 4.09 




w ^x 


3.41 




7.16 6.33" 



groups show a linear behavior with respect to Mendeleev table 
period of corresponding cation. Moreover, as shown in Table 
II and Fig. 2, bulk moduli show a less satisfactory agreement 
with experimental results ( maximum deviation within 20% 
in the case of II group, within 45% for the ji-PbFx). Also 
for bulk moduli of compounds with cations belonging to the 
II and IIB groups, an almost linear behavior results with re- 
spect to the Mendeleev's period. For IIB group no available 
experimental data supports the theoretical trend. 

Energy band structure, total and projected density of states 
(DOS)'" are shown in the region of the gap in Figs. 4, 5. Di- 
rect and indirect minimum band gaps are also clearly shown. 
In Table III these data are shown in comparison to previous 
results and experimental data; valence bandwidths are also 
reported in comparison with theoretical and experimental re- 
sults. In Table IV the main vertical transitions are reported 
for all the compounds. In the following lines, in order to ex- 
tract general trends for the electronic excitation properties in 
the stable cubic fluorides under study, we shall compare the 
results for the different fluorides considering first compounds 
with cations belonging to the II group, namely Ca, Sr and Ba, 
then compounds with metal belonging to IIB group, namely 
Cd and Hg. Finally for jS-P^A data will be analysed sepa- 
rately. All these compounds show a minimum direct band gap 



at the r point, except l5-PbF2. Concerning the smallest gap, 
the fundamental one, all the compounds herein treated are in- 
direct gap insulators apart from HgF2 which shows a direct 
fundamental band gap (Table III). Note that absolute values 
for the gaps show the so-called band-gap underestimate which 
can be resolved going beyond DFT-LDA, by using more ac- 
curate techniques for the excited states (e.g. GW ones, see 
Table III). This issue will not be addressed in the present 
work and will be the subject of future research. Moreover 
trends on electronic excitation energies as shown in the fol- 
lowing should not be affected by the above problem. '^ For 
Il-compounds, a decrease of the direct gap at F of 1.04 eV 
is similar to the 1.10 eV decrease after the experimental data 
(see also Fig. 4). On the other hand, a decrease of 0.54 eV 
for the X-F transitions for calcium and strontium fluorides re- 
sults which can be compared with a 1 .20 eV decrease from the 
experiments. Considering barium fluorite, the X-F transition 
is not the minimum, but that which occurs at (1/4,1/4,0)- 
F (see also Fig. 4c). However, considering the value for the 
X-F transition shown in Tab. IV, we confirm a smaller value 
of 0.31 eV with respect to the same transition for SrF2 (to 
compare with 0.6 eV experimental decrease). Considering the 
valence bandwidth, a 1.39 eV decrease going from Ca to Ba 
occurs, and that decrease can be compared to a 0.70 eV exper- 



Lattice Constant O 

Bulk Modulus © 

Experimental Lattice Constant v 

Experimental Bulk Modulus O 



Fundamental Gap D 

Valence Bandwidth O 

Experimental Fundamental Gap v 

Experimental Valence Bandwidth O 
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Figure 2: Lattice constant and bulk modulus for each fluorite 

under study. Results after present study Vs. available 

experimental data. 



Figure 3: Minimum energy gaps and valence bandwidths 

trends for the compounds under study. Also reported 

available experimental data. 



Table IV: Energy band gaps (eV) after present work. 



Direct gap 


CaF2 


SrF2 BaF2 CdF2 HgF2 p-PbF2 


L^L 


9.35 


10.59 


9.29 


7.62 


4.31 


4.62 


r^r 


7.71 


6.99 


6.67 


3.37 


0.35 


7.45 


x^x 


8.09 


8.03 


7.23 


8.18 


5.34 


4.09 


w ^w 


8.57 


8.44 


7.68 


8.17 


5.45 


5.10 


K^K 


8.50 


8.41 


7.54 


8.25 


5.34 


5.45 


Indirect gap 




x^r 


7.43 


6.89 


6.79 


3.04 






L^r 


8.00 


7.33 


6.73 


3.29 






w ^r 


7.54 


6.97 


6.90 


2.94 






K^r 


7.55 


6.98 


6.83 


3.14 






(1, i, o)^r 






6.58 








r ^x 












6.36 


w ^x 












3.95 


K^X 












3.98 


L-^X 












4.02 



imental one. 

For IIB-compounds (see Figs. 5a, 5b), the direct gaps at the 
r point show 3.02 eV difference going from CdF2 to HgF2. 
While the CdF2 presents an indirect fundamental gap between 
W — 7- F, HgF2 shows a direct fundamental band gap at F (see 
Fig. 5b). Moreover, for the IIB-compounds HgF2 presents a 
larger value of valence bandwidth of 0.59 eV. 

For l5-PbF2 the minimum direct gap occurs at X instead 
of at r as for the other fluorides. It shows an indirect fun- 
damental gap W ^ X of 3.41 eV (see Fig. 5c and Tab. III). 
Moreover, the largest value for the valence bandwidth occurs, 
with a slight overestimate of the experiment (13%). 

With respect to direct transitions, considering the data in 
Tab. IV, all Il-compounds show a decreasing trend at all main 
symmetry points, except for L. At that point SrF2 shows a 
1.24 eV larger band gap with respect to CaF2 and a 1.3 eV 
larger value with respect to BaF2. For IIB-compounds, the di- 
rect transitions show a larger difference (up to 3.31 eV) going 
from Cd to Hg. 

To complete the picture, the indirect gaps between the top 
of the valence band at high symmetry points to the bottom 
of the conduction bands (F for II- and IIB-compounds, X for 
li-PbF2) are shown in the second part of Tab. IV . 



IV. CONCLUSIONS 

The DFT-LDA electronic structures for ground-state and 
excited states for some cubic fluorides have been studied in 
detail. The electronic density of states (DOS) at the gap re- 
gion for all the compounds and their energy-band structure 
have been calculated and compared with existing experiments 
and previous theoretical results. The electronic energy-bands 
and transition energies are given and discussed. Within the 
same DFT-LDA scheme, general trends for the ground-state 
parameters and the DOS are also given. The above trends 
show good comparison with experimental data and theoreti- 
cal results. Relatively to electronic excitations, the conduction 
bands for Il-compounds are mostly dominated by the cations 




Figure 4: Energy band structures, total and projected density 

of states (DOS) are shown in the region of the band gap for 

the Il-compounds. 



d-orbitals, while for the IIB-compounds the tail of the DOS at 
the lowest conduction bands shows a largely s-type character 
The obtained DFT-LDA valence bandwidths agree with ex- 
perimental values within 30%. The present systematic DFT- 
LDA study is of particular interest for future researches on 
excited states and optical properties calculations of fluorides. 
We plan to carry out such calculations in next future by using 




those techniques particularly devoted to that issue. 



X W K 

(c) li-PbF2 



Figure 5: Energy band structures, total and projected density 

of states (DOS) are shown in the region of the band gap for 

the IIB -compounds and P-PbFj . 
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